New classes of spin chains from (SO( q )(N), Sp( q )(N)) Temperley-Lieb algebras: 
Data transmission and (q, N) parametrized entanglement entropies 
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A Temperley-Lieb algebra is extracted from the operator structure of a new class of N 2 x TV 2 
braid matrices presented and studied in previous papers and designated as SO( q )(N), Sf>( q )(N) for 
the g-deformed orthogonal and symplectic cases respectively. Spin chain Hamiltonians are derived 
from such braid matrices and the corresponding chains are studied. Time evolutions of the chains 
and the possibility of transition of data encoded in the parameters of mixed states from one end to 
the other are analyzed. The entanglement entropies S(q, N) of eigenstates of the crucial operator, 
namely the q-dependent JV 2 x iV 2 projector Po appearing in the corresponding Hamiltonian are 
obtained. Study of entanglements generated under the actions of SO(N), Sp(N) braid operators, 
unitarized with imaginary rapidities is presented as a perspective. 



I. INTRODUCTION 

Interesting developments in fault-tolerant quantum 
computation using the braiding of anyons have brought 
the theory of braid groups at the very core of topologi- 
cal quantum computing []_]. Moreover, the recent study 
by Kauffman and Lomonaco that the "Bell matrix", a 
specific braiding operator from the solution of the Yang- 
Baxter equation, is universal implies that in principle 
all quantum gates can be constructed from braiding op- 
erators together with single qubit gates [2j. In another 
very recent paper, the authors presented a new class of 
braiding operators from the Temperley-Lieb algebra [|[ 
that generalized the Bell matrix to multi-qubit systems, 
thus unifying the Hadamard and Bell matrices within the 
same framework Q. 

Here, we extract the Temperley-Lieb algebra from the 
operator structure of a new class of N 2 x A 2 braid matri- 
ces presented and studied in previous papers and desig- 
nated as 50(g) (iV), Sp( q )(N) for the g-deformed orthogo- 
nal and symplectic cases respectively Q . The connection 
between spin chains and Temperley-Lieb algebras is well- 
established Q. We derive the spin chain Hamiltonians 
from such braid matrices and study the corresponding 
chains. We then analyze the time evolutions of the chains 
and the possibility of transition of data encoded in the 
parameters of mixed states from one end to the other, 
following studies such as Ref. [_}. We further obtain the 



entanglement entropies S(q, N) of eigenstates of the cru- 
cial operator, namely the ^-dependent N 2 x TV 2 projector 
Po appearing in the corresponding Hamiltonian. Finally, 
the study of entanglements generated under the actions 
of SO(N), Sp(N) braid operators, unitarized with imag- 
inary rapidities, is presented as a perspective. 



II. BRAID MATRIX FORMALISM 



canom- 



In the context of exhaustive construction of 
cally factorized" forms [j| of braid matrices, new classes 
of solutions with remarkable and intriguing properties 
were obtained for SO q (N) and Sp q (N). To distinguish 
these special cases, they were respectively denoted by 
SO q (N), and Sp q (N). The aspects of that formalism Q 
necessary for our present study are summarized below. 

For SO q (N) and Sp q (N), the standard cases were first 
expressed [_| in the form 



R(6) 



/+(- 



U(0) 



-P+ 



/- M) 
/+(*) 



/o i-0) 
fo (0) 



Po, (1) 



where (P + ,P_,Po) form a complete basis of projectors, 
satisfying 



PiPn 



= ft • P 



P+ + P_ + Po = I. 



(2) 
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Since these are N 2 xN 2 gr-dependent matrices, so is R (6). 
All ^-dependence of R (9) is in the coefficients (/+ (9) , 
/_ (9), f (9)). As consequence of equation ^ 



R(9)R(- 



(n 2 xn 2 ) — In ® In 



I® I, 



(3) 
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where I is the (N x N) unit matrix. 

The factorization of the coefficients, along with ((2|), 
implies the "canonical factorization" of R (0) : 

R {9) = (/+ {-9) P + + /_ {-0) P_ + f {-6) P ) (4) 

x((/ + (e))~ l p + + (/_ (0))- l p_ + (/o (0))- 1 p o ). 

Like the projectors, the f's will also depend on q, the 
parameter of "(/-deformation" . 

For R {9)to be a braid matrix it must, by definition, 
satisfy the braid equation which provides matricial rep- 
resentation of the third Reidemeister move in the theory 
of classification of braids and knots. This means that 
defining (with "rapidities" {6,9')) 



R{9)®I ■ 
I®R{9) 



R12 {9) , 

i?23 (0) , 



we have 



B X =R 12 {9) i? 23 {9 + 9') R 12 {9') 
B 2 ee R 23 {9') R 12 {9 + 9') R 23 {9') . 



(5) 



The solutions for the coefficients must be found such that 
for a given set of explicitility defined N 2 x iV 2 dimensional 
projectors, one obtains 



B\ — B 2 



(6) 



In Ref. [5[, the coefficients for standard known solu- 
tions for SO q {N) and Sp q {N) were factorized as in (Q} 
and new solutions were obtained, which were studied in 
subsequent papers [9]. 

The new classes (denoted as SO q {N) and Sp q (N) re- 
spectively) correspond to 



/+ {0) = f- {O) = 1 
and hence (with new solutions for /o) 

R{9) = P+ + P_ + f °(~^ P 



(7) 



= J® J+( 



fo{6) 
fo{~ 



fo{0) 
I®I + lo{9)P . 



-l)Po 



(8) 
(9) 



In this paper, we analyze for the first time, the proper- 
ties of Pq that makes the solution of © possible. Later 
we consider u>{0) in that context. This turns out to be 
fruitful indeed. 



III. SO q {N), Sp q {N) AND TEMPERLEY-LIEB 
ALGEBRA 

We define 

{lu{9),lj{9'),lu{9 + 9')) = {lo,J,lo"). (10) 



Implementing © in ([5]), one obtains 

{Bi - B 2 ) = {u + u' + ujJ - uj"){{P ® I) - {I ® P )) 
+uu'u"{{P <g) I){I <g> P )(P (» I) 
-(/®P )(P ®/)(I®Po)). (11) 

We first define some notations: 

(1) i = N — i + I (when i = i). 

(2) {ij) as the (N x N) matrix with unity on row i and 
column j and all other elements zero. 



(3) The g-brackets 



[N±l} = 



q N±l _ q -N±l 

q-q- 1 



(12) 



We then have the projectors Pq as A^ 2 x A^ 2 matrices 



(1) For SO q {N), (iV = 3,4,5,...) 

JV 

{[N - 1] + 1)P = Q (pT ~ Pj) (ij) ® (y), (13) 
and 

(2) For5p g (iV), (iV = 2,4,6,...) 

Af 

([JV + 1] - 1)P = ?^~ W W((y) ® (*?))> (14) 
where 



= 1, (i < iV/2) 
= -1, (i > JV/2). 



(15) 



Note that we restrict q to be real, positive throughout 
so that one obtains real Pq. Also, note that the parame- 
ters p are A^-tuples: 

1. For SO{2n + l): 

1 3 11 1. 
p: {n- -,n -,0, -n+ -). (16) 



2. For SO{2n): 

p: {n- 1, n - 2, ...,1,0, 0,-1,..., -n+ 1). (17) 



3. For 5p(2n) : 

p : {n,n~ 1, 1, -1, -n). 



(18) 
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The projectors thus defined, can be shown to satisfy 
(P ® 7)(7 ® P )(P ® i) = fc~ 2 (P) ® J) (19) 
(7 ® P )(P ® /)(/ ® Po) - fc _2 (J ® P ), (20) 

where 

fc = ([JV - 1] + 1), for SO(iV) (21) 
fe = ([N + 1] — 1), for Sp(N). (22) 

This is the core of the Temperley-Lieb algebra to be de- 
veloped fully for spin chains, in a following section. 

At this stage, implementing (fl~9| and ([20|) in (fTT|) . one 
obtains 

(Pi -B 2 ) = (uj + uj' + ujuj' - uj" (23) 
+k~ Ww")((P ® /) - (/ ® Po))- 
Hence the braid equation is satisfied, if 

oj + J + ojuj' - uj" + k^ujuj'uj" = 0. (24) 

This non-linear functional equation was solved in our pre- 
vious paper @ . Denoting the special cases as SO and Sp 
henceforward, the solution is given by 

smn(77 + u) 



where 



and 



N-e _ -N+e 

(e" + e"" = ife = ? ^ + e (26) 

q-q 



e = 1 for SO (AT) 
e = -1 for SpfTV). 
We note the following points: 

1) Of the three projectors (P+, P_, pj), only one pro- 
jector, P , satisfies (fH?)) and (|2U)) . Instead of ©, if we 
start with 

R(0) = I ® 7 + AT±((9)P ± , 

then no solution is obtained, since P± do not satisfy the 
analogues of flU?]) and (j2"0")) . Hence, P is the crucial op- 
erator. 

2) An adequate solution of the non-linear functional 
equation in two variables (9, 6') is not evident to start 
with. But it does exist and accordingly, an explicit man- 
ageable solution (given by ([2"5j) and (|2"6")l ) is obtained. 

3) We can re- write (|2"6"|) as 

coshr, = h = ~([ATt1]±1) (27) 

(upper signs for SO(N) and lower signs for Sp(N)). We 
can also write 



sinhn = ±— V & 2 — 4 
' 2 



(28) 



for both cases SO(N) and Sp(N), and hence sinhr; can 
be chosen to be positive or negative, a point which we 
will revisit. 



IV. SPIN CHAIN HAMILTONIAN AND 
TEMPERLEY-LIEB ALGEBRA 

The chain Hamiltonian (and also higher order con- 
served quantities [5j) can be obtained directly as follows. 
We define 



R(O) = (—R(0)) e=o , 



(29) 



and with (R(Q))i t i + i acting on sites (I, I + 1), the Hamil- 
tonian for an open chain of length r can be written as 



77 = 7 ® ... ® 7 ® (P(0)) M+ i ® 7 « 



(30) 



i=i 



where 7 is the N x N unit matrix. 

For a closed chain with circular boundary conditions, 
there is an additional term with 



0R(Q))r,r+l, (r + lwl). 



(31) 



For our case 



2k 



(fc 2 - 4) 



T7^) P o (32) 



P(0) = -(2coth?7)P = =f( 



for the upper and lower signs in (|28p respectively. 
Hence for open r-chains, we have 



2k 



r-l 



H = T( , fc2 Z\y/2 I®-® I® (fl>)t,i+i ® J® ...<8i). 

(33) 

We define 

X; = J® ... ® I® (P )y+i ® 7 ® ... ® 7. (34) 

Then 



r-l 



P = -(2coth ?/ )(^^), 



(35) 



/=i 



where 77 is non-zero, and positive or negative according 
to the sign chosen in (l28l) . 

From (fT9" l) - (f22"l) with fc given by (|27]l . 



X? = Xi 
XiX m = X m Xi (\l - ml > 1) 



(36) 



Thus the chain Hamiltonian is obtained as a sum over 
generators of the Temperley-Lieb algebra, defined by 



Defining 



X[ = kX l = (e^ + e- r ')X' l , 



(37) 
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one obtains 

x l X l+i X 'l = x i 

X' 2 = (e" + e^)X[ 

X' l X! m =X' m X , l (|l-m|>l), (38) 

a standard form of Temperley-Lieb algebra. 
We note the following points: 

1. As we mentioned before, the link between spin 
chains and Temperley-Lieb algebras is a well- 
studied subject p. But here we have more than 
the defining relations (f36f or (|38|) . We have, for all 
N, explicit N 2 x N 2 matrix realizations of the gen- 
erators: (implementing ([T5 |) -([l8" ]) in (0D and (1371) 1. 
This, as will be displayed below, enables one to con- 
struct eigenstates and eigenvalues of chain Hamil- 
tonians for all N. 

2. The two signs in (|28|) will be seen to correspond 
to inversion of the sign of eigenvalues of H. They 
correspond to two different regimes. 

V. EIGENSTATES AND EIGENVALUES OF P 
AND ACTION OF H 

We start by presenting the action of the projector P 
on product states and then derive that of the Hamilto- 
nian ([55"]). The definitions (|12l) . (jT3"]l and the explicit 
particular cases of Appendix A, imply that the action 
of Po on general mixed states selects out specific linear 
combinations of states 

= \i) ® \N - i + 1) . 

In terms of 

^ = ([A- e ]+ e )P , (39) 
one obtains from (p ]) -([i"8 ]) . for Sd{N): 

N N N 

i^(E *a \a)) ® (J2 Vb \b)) = E 9 (p7_Pl) ^yi I") . (40) 

a=l 6=1 z=l 

and for Sp(N): 

N N N 

P^ xa |a» ® y » ^ = £ 9 (P7 " Pl) (e i ^)(£T2/7) |0) , 

a=l b=l i=l 

(41) 

with the e's defined in (fl"5"j). 

Explicit results of Appendix A can now be imple- 
mented as follows. 



where 

|*) = (<T 1/2 1 11) + |22) + q 1/2 |Tl». (43) 

The states (|1) , |2) , |l)) may be taken to correspond to 
spin projections (1,0,-1). This |\&) turns out to be an 
eigenstate of Pq with unity for eigenvalue: 

P |*) -I*) (fromP '|*) = (?- 1 + l + g)|*)), (44) 

(using fl39j) with N = 3, ([N - 1] + 1) = (g- 1 + 1 + g) ). 
From (|42]) 

P^(|lT),|22),|Tl)) = ( (Z - 1 / 2 ,l, (? i/2 ) | v|/) (45) 

and 

Po\ij) = {j^i). (46) 

Strictly analogous results hold for all SO(N) and Sp(N). 
This is already pointed out for the cases of Appendix A. 

B. SO (4) 



P^ = (q- 2 + 2 + q 2 )P . 

Then 

Po(xi\l) + x 2 \2) + %|2) + z T |l}) 
®( 2/1 |l)+t/ 2 |2)+y2|2>+y r |T» 

= OT^u/r + X2V2 + x iv-2 + qx-yi) I*) , (47) 

where 

|*) = q- 1 | LT> + 1 22) + |22> +g|Tl> (48) 

and 

P |*) = |*) . (49) 
As in (|45[) the non-zero results are obtained for only 
P^(|lT) , |22) , |22) , |Tl» = (q-\ 1, 1, q) |*> . (50) 

The states (|1) , |2) , |2) , |l)) correspond to spin pro- 
jections C§s 2 ' — 2 ' — §)• ^ ur notation here generalizes 
smoothly to any N. 

C. 5p(4) 



A. 50(3) 

Po(xi|l)+a; 2 |2) +x T |T)) 
(g)(i/ 1 |l) + 7/ 2 |2)+y r |T)) 
= (<T 1/2 x iy - + a; 2 y 2 + q 1/2 x Tyi ) , (42) 



P[ ) = (q- i + q- 2 + q 2 + q' l )Po. (51) 

Then 

iftii |l) + aj 2 |2) 1 2) + .t t |T)) (52) 
®(y 1 |l) + ?/2 |2) + y 2 -|2> + y r |T)) 
= (<2 , ~ 2 ziZ/r + q~ lx 2W2 - 3^52/2 - 9 2 a;Tyi) |*) , 
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where 



|*) = q- 2 1 11) + q- 1 |22) - q 1 22) - q 2 |ll) (53) 



and 



Non-zero actions: 



P |*) = |*) 



(54) 



P^(|ll) , |22) , |22) , |11» = (q- 2 ,q-\-q,-q 2 ) |*> 

(55) 

Let us now consider the action of H on an open chain of 
length r. 

Similar to ([3"9")l . Pq = fcP , we now define 



r-l 



X[ = kXi =^J® ...®I<g) (P^) M+1 ®I( 



/ (56) 



2=1 



with (/ = 1, r — 1), and 

k= ([AT T l]±l) 

9 - 9" 1 



±1, 



for SO{N) and Sp(N) respectively, where we consider 
real values of q. 

Using (f28|) . we can now write 



2k 



y/k 2 -A 
(sinhry) _1 P 



r—x ^ r — i 

1 sinh r\ 1 



1=1 
1 u' 



1=1 



(57) 
(58) 



Now 



P' \i x i 2 ...i r ) = ^ \hh— k-i) ((Pq) v) 

(59) 



and 



P(5 = 6(ii,ii+i)q^- p *hiei \im) (60) 



with the n-tuples defined by (|Tr?)l -(|18 p . and with e's as 
indicated by ([15]) for Sp(N) (each e = 1 for SO(N)). 

Let us concentrate on the explicit case of SO(3). Con- 
sider for example, the S*0(3) 4-chain with mixed product 
states 

|X) = (ai|l)+a 2 |2)+a T |T)) 
®(6 1 |1) + 6 2 |2) + ^|T)) 
®(ci|l)+c 2 |2)+c r |T» 
®(d 1 \l)+d 2 \2)+d T \T)) 
= \x) 1 (g) \x) 2 ® |x) 3 (g) |a;) 4 . 

Defining, as in (j43|) 



' 1 /2|iT) + |22)+g 1 / 2 |Tl) 



(61) 
(62) 



we now have 
P^ \x) = (a^q- 1 / 2 + a 2 b 2 + a^q 1 ' 2 ) |a:>3 |x) 4 
+ {b 1 c T q- 1 ' 2 + & 2 c 2 + bjCiq 1 ^ 2 ) \x) 1 |z) 4 
+(cid T q- 1/2 + c 2 cZ 2 + cydig 1 / 2 ) |x) 2 . 

(63) 

A generalization to a chain of length r is quite straight 
toward. In notations that are evident 

r-l 

H'(\ X ) 1 ... \ X ) r ) =J2fl - l*>J-l l*> l*>» + a - l»)r - 

(64) 

where 

/, = (a«4' +1 5 ? -V 2 + a « 6 (W) + ^0^+1)^1/2). (65) 

For 50(4), S , p(4) and so on, one can easily implement 
and (I59[) . with previous definitions. 



VI. TIME EVOLUTION OF SPIN CHAINS AND 
DATA TRANSMISSION 

A. Evolution in Time 

We are now in a position to start studying the evolu- 
tion in time t of a chain under the action of the operator 
e~ lHt . As often, we try to display some basic features by 
presenting results explicitly for a few restricted simple 
cases and indicating how to generalize them. 

Consider an 50(3) chain of spins, the projections for 
spin 1 being denoted as 



(|+),|0),|-)) = (|1),|0),|1». 



(66) 



From (|4"2"]) - (l4l)l) one sees then the iterative actions of P 
on the buildings blocks 



where 



I*)!*),!*)!*),!*)!*) \j) 



|*) = <T 1/2 1 11) + |22) + q 1/2 |Tl) 



(67) 



and take the values (1,0, 1), are the essential ingre- 
dients, along with the basic initial results: 

P^lT) , [22) , |Tl» = (q- 1 + 1 + <7)Po(|lT) , [22) , |Tl» 

= (,- 1 /2 ;1;(? l/2 )(g -l/ 2 | lI) 

+ |22)+ 9 1/2 |Tl» 
= ( g - 1/2 ,l,? 1/2 )|*) 



Po\ij)=0 {j^'i) 



and 



kP I*) = F |*> = (q- 1 + 1 + q) |*> = k |tf ) . (68) 
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From p8|). (|32J)-(|35J) and (EH]), for an open r-chain 



Next assuming j ^ i, we have 



H = \ ® ... ® I ® (Po)i,i+i ® I < 



(69) 



1=1 



where (corresponding to the sign of r\ chosen in (|28jl ) 

2 



A = T 



(70) 



which (from is real for q real, positive (which we 

assume to be the case). Note that for a closed chain the 
summation (|69[) would include an extra term I = r with 
(r+1) w 1. 
Defining 



if = \H' 



(71) 



we consider the series expansion (with I <£> I = Ig for 
SO(3)): 



-iff* -i\tH' 

e = e 

= I 9 + (-iXt)H' + ^(-iXt) 2 (H') 2 (72) 

+I(-at) 3 (JJ0 3 + - 

up to any chosen order in i. 
Suppose that the spin states 

|*)|i),|i)|*) s |i)|*) |j) 

correspond respectively to the sites 

(1,1 + 1,1 + 2), (1-1,1,1 + 1), (1-1,1,1 + 1,1 + 2). 

Defining 

H '(3) =Po®I + I®Pq 

they will always be implicitly assumed to correspond to 
the appropriate sites as the relevant parts of the total H' 
acting on the total chain. Thus HL^ \^>) \i) corresponds 

to the terms of H' acting on the sites (1,1 + 1,1 + 2) and 
so on. 

One obtains from ([32 ]) - ([15 ]) . 

^ 3) |*)|i) = fc|*)|i) + |i)|*) (73) 
H( 3) |i)|*) = fc|i)|*) + |*)|i). (74) 

Iterating, one obtains (see Appendix B) 

(Hl 3) n\*)\i))=A p \*)\i)+B p \i)\V) (75) 
(H' (3) y\i)m=A p \i)\*)+B p \*)\i), (76) 

where 

1 



.4 



p __((fc + l)f + (fc-l)f) 



b p = -((k + iy-(k-in 



(77) 
(78) 



H[ 4) (\i) W |j» = |*> \ij) + k \i) |*> |j) + |*) (79) 
or 

(H' w -k)Qi)\*)\j)) = \ij)\*) + \9)\ij). (80) 

Now, from (fTg ]) . ([74) and (|5U |) 

(i/ ( ' 4) -fc) 2 (|z)|*)| J )) = 2|z)|*)b-) (81) 

since Pq = 0, (J Hence 

(P ( ' 4) - fc) 2 "(K) |*) = 2(- 1 >(| ? ) |*) |i)) (82) 



(^(4) - kf n+1 {\i) |*) |j» = 2(«- 1 )(|ii) I*) + I*) |y> . 

(83) 

For |j) = |i) there are additional terms. One obtains 
#(4)(H) l*> |*» - l*> + * l») I*} \i) + I*) (84) 



H' (i) m |m)) = fc|*) |m) + |t) |*) |?) + |*> |*) (85) 



#(4)(|«) l*» = *|<?) I*) + K> I*) \^) + Q 5i l*> I*) - (86) 



where Si = (— ^,0, |) respectively for 



(1,2,1). 



(87) 



Hence, 



(fffo) a (|i>|*) |I)) = (fc 2 + 2)|z)|*)|I) 

+2fc(|*) | it) + |*)) 
+2q 5 < |*) |*) 



and again 



^ 4) (|*)|*))=2fc|*)|*)+X;g«'|i)l*) F) 

= 2fc|*)|*) + (<r 1/2 |l)|*) |T) 

+ |2>|*}|2}+g 1 / 2 |I)|*)|l}). (89) 



Using the set (|84]) - ([89|) one can now iterate. The way 
to proceed and the essential ingredients have been all pre- 
sented above. We will not write down the general result 
for (-ffm) n (|i) l^) |i))- In all the examples above there is 
one feature in common: The eigenstates | Vf) of Pg appear 
under the action of H' and move along the chain under it- 
erations. They move both forward and backward. There 
can be multiple |\&) depending on the initial state. The 
iterations above are to be implemented in (|72[) . Using 
systematically the results above one can start to study 
the evolution of an initial chain configuration. 
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B. Transmission of Data along a Chain 

For clarity and relative simplicity we start with an open 
6-chain of 50(3) spin states (f66]). The initial configura- 
tion is assumed to be (at t = 0) 



\X) 



(o) 



( Cl +c 2 |11» 11111) (90) 
Cl \X) l+ c 2 \X) 2 . (91) 



(We do not immediately normalize \X) for convenient 
generalization to more parameters, such as, that to start 
with (|99"]) . considered later). As will be shown below, 
time evolution under the action of eT %Ht will generate 
(at time t) a mutually orthogonal set of states including 



|llll) (di |ll) + d 2 |±T». 



(92) 



The other states at a finite non-zero t will be (apart from 
([90]) ) sequences 

{\X) X ,\X) 2 ,\V) |1111),|1)|*)|111), 

|11) |*) |11) , |111111) , |111221) , |111122)), (93) 

whose coefficients can be obtained (see Appendix C). 

With t increasing, the coefficients of the above set, (i.e. 
(f9"2"| and (f9"3"|) ) will continue to change. Otherwise, the 
set will be stable, no new basis states of the 6-chain will 
appear. This is a consequence of the specific properties of 
our H. It will be shown that it is sufficient to implement 
the series development 



-iHt 



i(\i)H' 



+ i(Ai) V) 4 - ^(Af)W + 0(t 6 ). (94) 

Evaluating finally, 

(H'r(ci\X) 1+ c 2 \X) 2 ) 

for p = (0, 1, 2, 3, 4, 5) one already obtains states of the 
type (|9"2")l (along with others orthogonal to it as given in 
(1931) ). Moreover (di, d 2 ) is obtained explicitly in terms of 
(ci,C2,A, t) where t is given (for any chosen origin) and 
from ((70)) (restricting the values of q for dcfiniteness, to 
q = 1 + S, 5 > 0) 



A 



V5 



(95) 



for 50(3) (i.e. for k = 3). The sign ambiguity in ([70)) 
corresponds to the two possible determinations of r/ (as 
explained in (|28p) corresponding to two possible regimes. 

Next, one can easily invert the relations and thus ex- 
tract (01,02) from (di,d 2 , X,t). Thus the initial mixed 
state at the left of the chain can be extracted by precise 
observation of the specific mixed state (|9"2")) at the right 
end of the chain at a finite time t. 



In this precise sense, we say that the initial state 
(ci J ll) + C2 I ll)) at the left has been transmitted to the 
right as (di |ll) + d2 where (01,02) can be recovered 
from (d 1? d 2 )- 

From the results of Appendix C one obtains 



di = C1X1 + c 2 x 2 
d 2 = CxVi + c 2 y 2 , 



(96) 
(97) 



where (x\,y\), (x 2 ,y 2 ) are given in the appendix. From 
(|179|) . (|181|) one sees (since A and k are known) that the 
coefficient of (t 3 ) in x 2 gives directly C2 from d\. One then 
easily extracts c\ also from the coefficients of powers of t 
in (di,d2). Thus our goal is attained. 

Apart from the development (|72[) in powers of t, we 
can also set 



q=l + S 



(98) 



and assuming 6 to be small, use a series development in 
powers of S to extract information more readily concern- 
ing the initial state from that at time t. Let us illustrate 
this, very briefly, using a simple example. 
Generalizing (|90|) to 



H'\ 



>(0) 



(a |Tl> + 6 122> + c |ll>) |1111) , (99) 

(q-^a + b + q^c)^) |1111) 

+ cg 1/2 |1) I*) |111) (100) 



((«- 



-5(a-c) + 6)|*)|llll) 



+c(l + -*)|l)|*)|lll) + 0(<5 a ). (101) 

Thus, in the corresponding generalizations of (JTO)), (|9T)) 
one can separate (a + c) and (a — c) and hence (a, c) by 
extracting coefficients of powers of 5. 

This can be more helpful for more elaborate initial 
states. A lesser number of powers of t will be needed 
to extract the initial parameters from the generalizations 
of (di,d2) above. 

The special significance of the point q = 1 will be em- 
phasized at the end of Sec. IVIIl in the context of entangle- 
ment entropy. Here we note that a supplementary series 
development about q = 1 can help in another context. 
The dependence of our model on the quantum deforma- 
tion parameter q is indeed a central feature. 

One may compare our results above with the study of 
"Quantum communication through an unmodulated spin 
chain" in 0. There one has only Pauli matrices and only 
two possible spin states. But the Hamiltonian couples 
all possible pair of sites and static magnetic fields are 
present. The action of e~ lHt is studied numerically. The 
specific structure of our Hamiltonian (not only as here, 
for N — 3 but also for N > 3 through straightforward 
generalization) make explicit computations feasible. 

To illustrate the above statements we consider, for 
50(3), a particularly simple initial configuration. 
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Suppose the chain C is given symbolically (with, 1, 2, 1, 
corresponding to spin projections (+1,0,-1) respec- 
tively) by 

Co = (...llll (p) T (p+1) m...) (102) 

with all sites up to, say p in state 1 and then all sites in 
state 1. 

Gathering together the definitions and notations (j5"6")) - 
(|59| in the compact notation 

(-iHt) = C[(-P (12) ® I ® I ® ■•■) + (J" ® ^0(23) ® ^ ® ••■) 

+(/ ® / ® P 0(3 4) ® -.) + •••] (103) 
= J2C(H p , P+ i). (104) 

p=0 

To start with, only H p ^ p+ i will have a non-zero action 
on C. 

HC = ...lll(* PiP+1 )TTT..., (105) 

where ^p.p+i (given by P5|) ) is 

(g- 1 / 2 |lT) + |22)+ g 1 / 2 |Tl)). (106) 

In the action of e~ lHt , the ^ states then spread out as 
follows 

{-iHtf : 

H% lP+ i -J- (*p-i, p , *p, P +i, * p +i,p+2) (107) 

(-iffi) 3 : 

H^ p ,p — > (^j,_2,p-i,^p-i,p,^p,p4-i) 
-ff^p.p+i - ^ (^p-i.p) ^p.p+i) ^p+i.p+a) 

#*p+l,P+2 (*p,p+l)^p+l,p+2)*p+2,p+3)- (108) 

One already sees, schematically, how the "J states are 
generated, move forward and backward, crossover and 
acquire coefficients corresponding to diferent terms 

— (-iHt) n (n = 1,2,3). 

Already the multiplicities counting the contributions 
from different ri's (with appropiate coefficients for each 
order) are for 

(*p-2,p-l, *p-l,p, ^p.p+l, ^p+\,p+2, ^fp+2,p+3) 

(1,3,5,3,1). (109) 

After, say r steps (i.e. up to order t r ), the above sequence 
is 

(..., 2r - 5, 2r - 3, 2r - 1, 2r - 3,2r - 5, ...). 

An initial state less simple can give a much more complex 
pattern. For initial C 

C+_ : (11111.. .Ill) 
C_+ : (TlTlT...lTl) 
C 00 : (22 2 22. ..222), 



even the term {—iHt) in the expansion of e~ lHt generates 
^ states for each pair of sites (p,p + 1). 

These states ^> are the basic building blocks of our for- 
malism. They will be seen to be entangled states, and 
their (q, A r )-dependent entropy will be studied in the fol- 
lowing section. 

VII. [q, iV)-DEPENDENT ENTANGLEMENT 
ENTROPY OF EIGENSTATES OF P 

Acting on the pure product states |m), (i = N — i + 1) 

for both SO(N) and Sp(N) the projector Pq creates its 
eigenstate 

N 

P o|«)=E« P o)(^),(7I))l^)«l VI/ )' ( U °) 
j'=i 

the matrix elements of Pq in (|110[) , being defined as in 

(HI - ©. 

Does Po thus generate entanglement? We give an affir- 
mative answer below, evaluate the entanglement entropy 
to quantify it and analyze the (<?, N)- dependence, q being 
the parameter of quantum deformation. We formulate 
the g-dependence, separately for different values of N. 

A. Sd(3) 

As usual, we start with 50(3) (see (|32 ]) -([l6" l) ) and 
(|135|) - (|138|) ) and as essential step, normalize \ty) as be- 
low, denoting it now by \^)/ n ) ■ Define 

l*)(„) = (Q- 1 + l + q)- 1/2 (q- 1/2 |+-) + |00) + q^ \-+)) 

= ( 9 - 1 + l + (Z )- 1/2 ( (Z - 1/2 |lT) + |22) + 9 1 / 2 |Tl)) (111) 

That this is an entanglement state, is evident immedi- 
ately. Attempting to re-express it as a product state 

(c+ |+> + Co |0) + c_ \-))(d+ |+) + d |0> + d. |-)) 

one runs directly into contradictory constraints on (cj, dj) 
coefficients. 

To quantify the entanglement one first notes that the 
eigenstate |^)(„) satisfying 

fi> !*>(„) = I*>(n) ( 112 ) 

is already Schmidt decomposed as 

l*>(») = $>«!*>• ( 113 ) 

i 

Hence, without passing via the density matrix (and with- 
out using a log 2 basis) one obtains the von- Neumann 
entropy [lfl[ll| as 

S=-EM 2ln M 2 ' i=(+ s 0,-)- (H4) 
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From pTTI) . (|TT3l) and ATM]) , one obtains 
5(g) = S(q~ l ) 

= i ^" 1+i+ ^- (g (9 r + CV 119 - 



(115) 



From g = 1 one obtains the maximum entropy as 

5(1) =5(max) =ln3. (116) 
To first order in e > 0, 



5(1 T e) = In 3 ± -eln(l T e) < 5(1). 
o 



(117) 



Consistently with (fTT5|) (namely 5(g) = 5(g -1 )) 

5(g) -4- as q -> oo or g -> 0. (118) 

After displaying the q- dependence for N — 3 we explore 
below also the TV-dependence. As a hrst step we move 
up from N = 3 to N = 4. 

B. SO (4) 



From ([47 |I -([S0 | and (|T40l) - (fT42| we define now the nor- 
malized eigenstate of Pq as 

!*)(„) = (q- 2 +2+q 2 r 1/2 (q- 1 |lT) + |22) + |22)+g |Tl». 

(H9) 

The corresponding entropy is obtained as 
S(q) = Siq- 1 ) = 2 Hq- 1 + q) - 2 ^_|j In q. (120) 

Now, again for q = 1, 

5(max) = 5(1) = In 4 (121) 
5(1 Te) =ln4±4eln(lTe) < 5(1). (122) 

Once again ()118p 

S{q) ->• as g ->• oo or g ->• 0. (123) 

Hence as iV increases from 3 to 4, 5(max) moves up 
from In 3 to In 4 and falls a bit more steeply, but again 
symmetrically in (g,g _1 ) to vanishing asymptotic values 
as q — ► oo and q — > 

C. Sp(4) 

For iV = (4, 6, 8, ...), i.e., for each such even N, one has 
the projector P for Sp(N) as well as for SO(N). Though 
this paper is mostly devoted to a detailed study of 50(3), 
after showing how the entropy depends on N for SO(N) 
by presenting the results for 50(4), we also present 
briefly the results for Sp(4) to display both the analo- 
gies and the differences in this respect between 50(4) 



1 1 

1 2 
1 

in 8 
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FIG. 1. (q, iV)-dependent entanglement entropy. 

and Sp(A). The relevant generalization for N = (6, 8, ...) 
is straight forward. 

For 5p(4), starting with (|5i" l) -([54" P and ([Tg3 ]l -([Ti5 ]l one 
obtains (as compared to (|119p ) 

l*) (B) = (g- 4 + g- 2 + g 2 + g 4 )- 1/2 (g- 2 |lT) 

+g- 1 |22)-g|22)-g 2 |Tl)), (124) 

the corresponding entanglement entropy is 

5(g) = 5(g- 1 ) = ln(g- 4 + g- 2 + g 2 + g 4 ) 



4g 4 + 2g 2 - 2g~ 



4g" 



g^ 4 + g- 

Thus again, as for 50(4), 



g 2 + g 4 



■ In g. 



(125) 



5(1) = 5 (max) =ln4 

5(1 Te) = ln4±10eln(lTe) 
5(g) — > as g — > oo or g — > 0. 



(126) 
(127) 
(128) 



The slope, starting from 5(1) towards the asymptotic 
zero values is steeper as compared to the 50(4) case, as 
shown in Fig. [T] 

One can show (starting with Pq defined in (TT3l) . (TT4)) ) 
that for each N, 



S(q,N) = S(q-\N) 
5 max (g,A) = 5(1, N) = In N 
5(g, N) — > as g — > oo or g 



0. 



(129) 



We do not present the explicit (straight forward) compu- 
tations for general cases. But one aspect is worth point- 
ing out. 

We emphasized in related previous papers [5J, |9| that 
for our special solutions (SO(N), Sp(N)) the projector 
Pq and the braid matrix R remain non -trivial for q = 1 . 
This is a remarkable feature (absent in standard solu- 
tions for SU(N) q , SO(N) q , Sp(N) q ) as was emphasized 
in Section 3 of reference Si and Section 2 of reference M . 
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Now we have found that not only SO(N), Sp(N) re- 
main non trivial for q = 1, the associated entanglement 
entropy (for eigenstates of Po) is maximal (In TV) for 
9=1. Thus the striking non-triviality at q = 1 acquires 
further significance. 

An adequate study of correlations in presence of multi- 
ple states |\&) (as already starting to appear in (|88|l) will 
not be undertaken in this paper, this aspect remains to 
be explored. 

Wc have shown above that the operators Po, acting on 
any either annihilates it (for j ^ i) or generates (for 
j — i) entangled states and we have quantified the 
corresponding entanglement for all (q, N) by computing 
the entanglement entropy. 

VIII. CONCLUSIONS AND PERSPECTIVES 



tary by implementing imaginary rapidities (i0,i9') in R 
matrices of ([5]). Here again (from ([5]) and ©) 

R m = p + + p. + i^lp 

Jo (w) 

= I®I + w(i6)P (133) 
can be directly verified to satisfy unitarity, i.e. 

(R q (i6^R q (i0))=I N ^I N . (134) 

Now one can try to formulate explicitility (g, N)- 
parametrized entanglement quantifiers of the superpo- 
sitions of 3-qubit states generated by the action of the 
braid operator (see (0, ((U) on such product states, as 
on the l.h.s. of (|40p . (1TT1) generalized to triple products 
One can also examine possible teleportation protocols as- 
sociated to our class of unitary matrices (see Ref. 0]). 



Starting with the projectors (P+,P_,Po) for SO q (N) 
and Sp q (N) braid matrices and then keeping only P it 
was shown how Po can generate a Temperley-Lieb alge- 
bra and how this property leads to a remarkable special 
class of braid matrices (denoted as SO q (N), Sp q (N)) and 
related spin chains. 

Then we have explored certain aspects of such spin 
chains, using mostly SO q (3) examples of chains with free 
ends to display some particularly interesting properties. 

Time evolution of such chains was studied by evaluat- 
ing the actions of successive terms (—iHt) p in the series 
development of e~ lHt , H being the spin chain Hamil- 
tonian. In particular, we studied in what form the data 
encoded in parameters of mixed states at one end of the 
chain can be decoded by observing mixed states reaching 
(as t increases) the other end of the chain. Most of the 
relevant computations has been collected together in the 
Appendices. 

Finally, we have obtained the entanglement entropies 
S(q,N) of the eigenstates of Po. In particular we ob- 
tained (q, A^-dependence as 

S(q,N) = S(q-\N) (130) 
SmUliN) = S(1,N) = In N (131) 
S(q,N) -> as q^ oo or q^ 0. (132) 

We pointed out before in Sec. IIVI that the two possible 
sign determinations of the essential parameter rj corre- 
spond to two different regimes for the energy eigenvalues 
of the chain Hamiltonian. One may compare and con- 
trast such a feature with the well-known corresponding 
ones of the 6- vertex models (see for example, |12|V 

Certain aspects of our classes of spin chains remain to 
be studied, as pointed out in Sec. IIVI 

Another rich perspective is the exploration of various 
aspects of the braid matrices we started with (Sees. [II] 
and IIIII) before extracting from them the chain Hamilto- 
nian (Sec. |IV|) . 

In previous papers (T3 , [Tij we studied parametrized en- 
tanglements generated by braid operators rendered uni- 



IX. APPENDIX A 
Explicit P (N = 3, 4) 

Many basic results of Sec. [V] can be read off easily 
from the matrices Po presented below. The matrices (ij) 
are defined above (fT2|) as are (ij). The projectors Po 
are defined by (|12l) - (fT5|) . Their contents for the simplest 
cases are displayed below. 

i) SO(3): (N = 3; 1= 3, 2 = 2) 

(q~ l + 1 + q)P Q = P^ 

= g -1 (ll) ® (TT) 

+q- 1/2 (12)<E> (12) + (11) ® (11) 
+q- 1/2 (2l) <g> (21) + (22) ® (22) 
+q 1/2 (2T)<Z> (21) + (Tl) ® (11) 
+ 9 1/2 (12) ® (12) + q(Tl) ® (11) 

(135) 
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Defining the base states 

\ij) = \i) ® |j) , 
the single eigenstate of Po with non-zero eigenvalue is 

|*) EE q- 1 ' 2 | IT) + |22) + q 1 ' 2 | Tl) , (137) 
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Po I*) = I*) • (138) 
All 6 states with j ^ i are annihilated by Po- Also, 
P (<7 1/2 |11) - |22» = P (|22) - q- 1 ' 2 |Tl» = 0. (139) 

Corresponding patterns arise for all iV. They will not be 
explored in any detail. 

We present briefly the cases N = 4. 

ii) SO(4): ( 1 = 4, 2 = 3) 

(< 7 - 2 + 2 + (? 2 )P = P ( 5 

= cT 2 (ll) ® (TT) + g _1 (12) ® (12) 
+g- 1 (12)® (T2) + (TT)® (11) 
+9^ (21)® (21) + (22)® (22) 
+(22) ® (22) + g(2T) ® (21) 
+g~ 1 (2T)® (2T) + (22)®(22) 
+(22)® (22) + g (2T)® (21) 
+(11) ® (11) + q{12) ® (12) 
+g(T2)® (12) + g 2 (TT)® (11). 



Defining 

|*) = -T 1 | IT) + |22) + |22> + g|Tl) , 

Po |*> = |*> • 
Hi) Sp{A): (1=4,2 = 3) 



(<f 



q- 2 +q 2 + q A )Po = 



q- i {ll) ® (TT) + 9 ~ 3 (12) ® (12) 

- g - 1 (l2)®(T2)-(lT)®(Tl) 
+g- 3 (21) ® (2T) + <T 2 (22) ® (22) 
-(22)® (22) - g(2T) ® (21) 
-^(fl) ® (2T) - (22) ® (22) 
+g 2 (22)® (22) + g 3 (2T)® (21) 

-(Tl)® (IT) - g (T2) ® (12) 
+ g 3 (T2)® (i2) + g 4 (n)® (ii). 



Defining 



P |*) = |*). 



(140) 

(141) 
(142) 



(143) 



|ll) +q^ 1 22) - q\22) - q 2 |ll) , (144) 



(145) 



For Sp(N) the blocks with negative signs are anti- 
diagonally aligned. 



X. APPENDIX B 



Iterative Action of H 



As explained in Sec. IVI| in studying 50(3) chains, it 
is useful to have ready results for 



and 



tf&(|i>|*> li))- 



Here we collect the results indicating the derivations. 
We consider below HLy HL, as defined below ((72)) in 
"notation" . Neccesary multiplicative factors can be eas- 
ily supplied. We start with results ff73|) . ((74)) . Using them 
one obtains 



ff ( ' 3) (|*)|i)) = (P '®7 + 7®P ')(|*)|i)) 

= k\*)\i) + \i)m, 

where k = (q^ 1 + 1 + q), i = (1, 2, 1) also 



(146) 
(147) 



Iterating 



(h( 3) ) p (I*)K)) = ^I*)I<) + 5pK)I*) ( 148 ) 

(ff ( ' 3) r(|i) |*)) = A p |i) I*) + B p I*) |i) , (149) 



where 



• 1 • 



((fc + l) p + (fc-l) p ) (150) 
((fc + lf- (fc- l) p ). (151) 



Now 



H{ A) {\i) |¥) |j)) = (i* ® / ® I + J ® P^ ® / (152) 

+/®/®P^)(|i)|*) |j)) 
= |i)|i)|*) + |*)|£)|i) (153) 

+*|i>l*> U) 



or 



(7i ( ' 4) - fc)(K> |*) |i)) = \i) |i) I*) + |*) |i) |i) . (154) 



For j ^ i (when i =^ j), Pg = and from (|146|) . 

dUTji, (irki) 



(PT ( ' 4) -fc) 2 (K)|*) |j» = 2|t)|*) |i) 



(155) 



Thus 



(/f ( ' 4) -fc) 2 «(|i)|*)|i)) = 2"- 1 |i)|*)|i) (156) 
(#(4) - k) 2n+1 {\i) |*) |i)) = 2- 1 (/P ( ' 4) - k) \i) |tt) li) 

= 2- 1 (| J ) |j)|*) + |*)|t) |i)). 
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These results can be implemented directly by writing 

e -iHt _ e ~ikt e ~i(H-k)t 

and using the series development of the last factor. For 
j = i there are extra terms as follows 

(flfo - k)(\i) |*) |5» = \iT) |*) + |*) . (157) 

Hence 

(PT ( ' 4) - fc) 2 (K) |*) \i)) = 2 \i) |*) \i) + 2 |*> |*> . (158) 
But now 

H' {i) |*> |*) = 2 |*) |*) + (/ ® P ' ® 7) |*) |*) 

= 2|*>|*} + (<r 1 /2|l)|vlr) |T) 

+ |2)|vI/)|2)+ 9 1 /2| T )| v p ) | 1))i (159) 

P ( ' 4) (<T 1/2 |1) |*) |T) + |2) |*) |2) + gV3 | T ) |1» 
-2|vI/)|vI/) + fc( g - 1 /2| 1) | vl/) |i) 

+ |2)|vI/)|2)+g 1 / 2 |T)|*)|l)). (160) 
Combining (jl 59[) . (|1 60[) one can now iterate. 



XI. APPENDIX C 
Explicit results for a 6-chain 

We present below the iterated action of H' (up to fifth 
order, namely (H'f) on the free 6-chain states (for our 
SO(3) case) 

13^ = 1111111) (161) 
\x) 2 = |lTllll). (162) 

They will be implemented in Sec. I VII to study, explicitly 
for a simple case, the time evolution of our class of spin 
chains and possible data transmission with such evolu- 
tions. 

Here the relevant H' is (with Pq defined in Sec. Iv| 

5 

H' {6) = ^ I ® I ® ... (Pq) m +i ® - ® P (163) 

The actions of Pq on 50(3) states are defined in Sec. 
IVl and iterative actions are presented, for H' acting on 
50(3) states in Sec. IVII For the restricted case relevant 
here only one needs, for sub-chains of P( 6 ) above 

H' (3) =P^(E>I + I®P^ (164) 
P (4) =P{ ) (E)I®I + I®P{ ) ®I + I®I®P{ ) , (165) 



acting respectively on 

flfo(|*>|l),|l>|*» (166) 

= ((fc|*)|l) + |l)|*)),(|*)|l)+fc|l)|*))) 
H' w (\l) |1)) = |*) |11) + k |1) |*) |1) + |11) |*) . (167) 

Here we have used the basic definitions and results 

(Sg)-®. 

Using all these results systematically one obtains the 
following results in a straightforward fashion, arranging 
terms in the order shown below. 

H'\x) x = g 1 / 2 |^> Jllll) (168) 
H' \x) 2 = q- l H' \x) 1 + q 1 ' 2 |1) |*> |111) , (169) 

where = q- 1 ' 2 |lT) + |22) + q 1 ' 2 |Tl). 

(H'f \x), = q l '\k |*) |1111) + |1) |*) |111)) (170) 
(H'f \x) 2 = q~\H'f |x) x + ^(l*) |1111) (171) 

where k = ((p 1 + 1 + q). 

(P') 3 ^)i=9 1/2 ((fc 2 + l)l*> inn) (172) 
+2fc|l)|*)|lll) + |ll)|*)|ll>) 

(H'f \x) 2 = q-\H'f \x) x + q 1/2 (2k |*) |1111> 

+ (fc 2 + 2)|l)|*)|lll)+2fc|ll)|*)|ll) 
+ |H1)|*>|1)). (173) 



(H'f \x) 1 = q^ 2 ((k 3 + 3k) I*) |1111) (174) 

+(3fc 2 + 2)|l)|*)|lll) 

+3fc|ll)|*)|ll) + |lll)|*)|l)) 
{Hr\x) 2 =q-\H'f\x) 1 (175) 

+(7 1/2 ((3fc 2 + 2)|*)|llll) 

+(fc 3 + 6fc)|l)|*)|lll) 

+(3fc 2 + 3)|ll)|*)|ll) 

+4fc|iii)|*)|i) + |mi)|#)). 



[H'f \x) 1 = q 1 / 2 ((fc 4 + 6fc 2 + 2) |*) |1111) (176) 
+(4fc 3 + 8fc)|l)|*)|lll) 
+(6fc 2 + 3)|ll)|*)|ll) 
+4fc|lll) |*) |1) + |1111) |*)) 



(H'f\x) 2 = q- 1 (H'f\x) 1 (177) 
+g 1 / 2 ((4fc 3 + 8fc)|*)|llll) 
+(fc 4 + 12fc 2 + 5)|l)|*)|lll) 
+(4fc 3 + 13fc)|ll)|*)|ll) 
+ (7fc 2 +4) |111) |*> |1> + 5fc |1111> |*)). 
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From the preceding results the coefficients (up to 0(t 5 )) _ _L(;^) 3 _|_ — (\t) 4 (4k + o + q^ 1 ) (181) 

of the states 1 1111) (|Tl),|22), 1 11)) are obtained as 3! 4! 

given below __L (A t)5 (7fc 2 + k{5q + 4g -i } + g) + ^ 6) 

-iAiif' ijs^v = ... + Mill) (ari 111) (178) 1 ' . i 

e- lXtH ' \X) 2 = ... + |1111) (x 2 |Tl) (179) z 2 = -l i (\t) 5 q 1 / 2 (5k + q~ 2 ) + 0(t 6 ). 



and 



+lfe|ll)+a&|22», 



x x - i(At) 4 - ^(At) 5 (4fc + q) + 0(t e ) (180) 

?A = -^(At) 5 + o(i 6 ) 
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